A periodic orbit basis for the Quantum Baker Map 
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A set of quantum states, dynamically related to the classical periodic orbits of a chaotic map, is 
used as a basis in which the description of the eigenstates of its quantum version is greatly simplified. 
This set can be improved with the inclusion of short time propagation along the stable and unstable 
manifolds of the periodic orbits resulting in a construction similar to the scar functions of Vergini 
[Jour. Phys. A 33, 4709 (2000)]. The average participation ratio is used to quantify the quality of 
the basis. 
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I. INTRODUCTION 



The question of the structure of the eigenfunctions of 
chaotic systems is intimately related to the construction 
of the bases in which they can be expanded. In this re- 
spect a good basis is one in which most eigenstates are 
given in terms of a limited number of significant coef- 
ficients: a good measure being i.e. the average partic- 
ipation ratio in the given basis. In the generic case of 
a basis unrelated to the Hamiltonian (or the map) the 
states are mostly random and the average participation 
ratio ((PR)) takes the random matrix value D/2 (for a 
Hilbert space of dimension D). At the other trivial ex- 
treme the eigen basis gives a (PR) of unity. It is then 
clear that a "good" basis has to incorporate some dy- 
namical elements from the system and at the same time 
be sufficiently simple so as to be effectively constructed 
without resorting to a full diagonalization. 
In the case of the quantum baker's map (QBM) Lak- 
shminarayan found that the eigenfunctions had a sim- 
ple structure (and significantly small participation ratio) 
when looked upon in the Hadamard basis [1], thus ex- 
ploiting a very special property of the QBM. This line 
of research was followed in [2], where it was found that 
the eigenfunctions of a large family of quantizations of 
the QBM could be described in terms of a very simple 
map, the essential baker, which for special values of the 
Hilbert space dimension becomes the Walsh quantized 
baker and can be explicitly constructed. These features 
are very special and intimately related to the binary sym- 
bolic dynamics of the map and are not easily generalized 
to other systems. 

A different approach with a potentially more general ap- 
plicability, is based on the construction of states that 
"live" on the unstable periodic orbits of the system. 
When used as a basis these states realize in quantum me- 



chanics the ideal of Poincare in the sense that "they (the 
periodic orbits) are the only breach through we might 
try to penetrate into a stronghold hitherto reputed unas- 
sailable" [3]. The fact that some eigenstates of chaotic 
systems show "scars" of periodic orbits was established 
long ago by Heller [4], in counterpart to earlier works in 
which the assumption was a uniform distribution on the 
energy shell according to the microcanonical ensemble 
[5, 6]. The scarring phenomena was studied in several 
chaotic systems, in which linear and non-linear theories 
were developed 0, i, i, ED, El El . 
The construction of scar functions in the QBM is based 
on early studies in the stadium billiard [TO, E2, E3 • The 
QBM can be thought of as a pedagogical system to apply 
these techniques since it has symbolic dynamics, stable 
(unstable) manifold parallel to q (p) direction, a finite 
spectrum and the same small- valued Lyapunov exponent 
in the entire phase space. 

In this paper we provide a recipe to construct a set of scar 
functions as an accurate basis to describe the QBM. This 
basis has the propagation time of the quantum propaga- 
tor as parameter. When this time is of the order of the 
Heisenberg time, the basis converges to the eigen base of 
the map. For short times, of the order of the Ehrenfest 
time, the basis describes the spectrum of the QBM bet- 
ter than other known bases p], Q • 

In section II, we briefly introduce the classical and quan- 
tum version of the baker map. We construct the periodic 
orbit modes and scar functions based on the evolution of 
the coherent states under the QBM. We give the rule with 
which we choose a basis to describe the spectrum of the 
QBM for any even dimensional Hilbert space in section 
III. Then, in section IV, we numerically test the basis, 
computing the average participation ratio as a function of 
the propagation time. In section V we propose a method 
to approximate the scar functions by homoclinic peri- 
odic orbit modes avoiding evolution in time, and finally, 
we state the conclusions. 
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II. CLASSICAL AND QUANTUM EVOLUTION 

A. The baker's transformation 

In this section we review some properties and nota- 
tion of the classical baker's map that we will use in the 
quantum states construction. The baker's map B [15[ is 
defined in the unit square phase space [0, 1)) as 



q' = 2q- [2q\ 
(p+[2q\) 



P 



(1) 



where [x\ is the integer part of x. This map is area- 
preserving, uniformly hyperbolic with Lyapunov expo- 
nent (A = In 2), and has stable foliation {q = est} and 
unstable foliation {p = est}. 

The baker map has a simple action upon symbols in 
the binary expansion of the coordinates 



(p\q) 



. 6-i • eoei 



(pf\qf) 



(2) 



where q = e l 2~^ and p = 

The map has two symmetries: 

• Parity (R): represented with the exchanges 
q —> 1 — q and p — > 1 — p together with the bitwise 
logical NOT upon symbols (0^1). 

• Time reversal (T): represented with the exchange 
p ^ q together with reversing the direction of the 
symbolic flow. 

The periodic orbits of the baker map of period L can 
be represented by binary strings v of length L. We de- 
note the different trajectory points on a periodic orbit by 
(qj,Pj) for j = 0, ••• ,L-1 with (q L ,PL) = (qo,Po)- The 
coordinates of the first trajectory point on the periodic 
orbit can be obtained explicitly in terms of the binary 
string 



v 



q — vvv . . . = 



Po = v^v^v^ . 



2 L - 1 
i/t 



2 L - 1 



(3) 
(4) 



where v is the integer value of the string v which repre- 
sents a binary number, and is the string formed by all 
L bits of v in reverse order. The other trajectory points 
can be easily calculated by iterations of the map or by 
cyclic shifts of v. 



B. The Quantum Baker Map 

The quantization of the map is performed in an even 
D-dimensional Hilbert space with D = l/(27rfi). The 



QBM is defined in terms of the discrete Fourier transform 
with antisymmetric boundary conditions as [l6|, [l?], EH 



B = G 



D 



G 



D/2 

G 



(j\G D \k) 



exp 



D/2 



(5) 
(6) 



The quantum baker map has the same symmetries as 
its classical counterpart 



B,R 







(7) 

(8) 



with parity represented by R = — G 2 and time reversal by 
T = KG, where K is the complex conjugation operator. 

The QBM spectrum is characterized by D eigenphases 
and eigenstates B\ipj) = e^'l^j), with definite R sym- 
metry (R\ipj) = and satisfying the time reversal 
requirement G\ipi) = \ipi)* . 



III. BASIS CONSTRUCTION 

A. Periodic Orbit Modes 

The first step in our construction is the definition of 
the periodic orbit modes (POM), a superposition of co- 
herent states centered on the periodic points of an or- 
bit. Similar constructions have been employed before 
B B EH El EH- Our definition is equivalent to the 
discrete-time version of the tube functions defined by 
Vergini and Carlo in [TO, [H| for continuous Hamiltonian 
flows. 

The coherent state on the D-dimensional Hilbert space 
on the torus with anti-periodic boundary conditions cen- 
tered on (q,p) ([20]) are represented in coordinate basis 
\j) as 



(ej ;+m — q/2)p— mm 



(j\q,p)=K J2 e-* D ^ +m -tfj 2 * D 

m= — oo 

(9) 

where ej — (j + 1/2) /D and K is a normalization factor 
which converges to (2/D) 1//4 for D ^> 1. The phase has 
been chosen in such a way that parity and time reversal 
operators act on them as 



R\q,p) = \l-q,l 
T\q,p) = \p,q) 



■P) 



(10) 

(11) 



without additional phases. 

We now consider the collection of coherent states on 
the periodic points \qi,Pi), i = 0, • • • , L — 1 of a given 
primitive orbit labeled by the binary string v. For chaotic 
systems the points are isolated and therefore in the semi- 
classical limit D —> oo these states are approximately 
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Figure 1: (color online) Top: Husimi representation of the 
periodic orbit modes over 

00! I^ob?"") ( left ) and 

over 00101 

|3>ob?oi + ) f° r D = 120 in the unit square phase space q,p G 
[0, 1). Bottom: Square modulus of the product of these states 
Mj\$ooT~)\ 2 (left) and \WA*lb*oi~)\ 2 with eigenstates |^> 



yj i^ooi 
with j — 1 



3 1^00101 

D ordered by increasing eigenphases. The red 



dashed line represents the central phase A 001 and A 00101 . 



orthogonal. In the same semiclassical limit they satisfy 
the approximate conditions 

(Qj+uP3+i\Qj,Pj) - S j+U ( 12 ) 

Bj j+i = (Qj+uPj+i\B\Qj,Pj) - 1 — = 

Vcosn A 

where qL = qo, PL = Po] A is the Lyapunov exponent and 
where the phase Sj acquired by the coherent state in one 
step of the map (with the present choice of phases for the 
coherent states) is 



Sj — SqjiPj 



1 



(13) 



The L x L matrix B 3 & in Eq. [12] is cyclic in the semi- 
classical limit and therefore it can be diagonalized by a 
discrete Fourier transform. The eigenvalues are given by 



.i27rA* 



%\B\<&) 



\J cosh A 



where 



At 



DS U + k 



(14) 



(15) 



The phase of the eigenvalues involves the classical action 
of the orbit S u = ^2f=o Sj, and k is a Bohr-Sommerfeld 
like parameter which can be chosen from fc = 0, . . . , L— 1. 



Each periodic orbit then contributes L complex eigenval- 
ues whose phases are equally spaced and shifted from the 
origin by DS U /L. 

The eigenfunctions are the periodic orbit modes. They 
are given explicitly by 



L-l 



\4>i) 



exp 



.2ir(DS v + k)j 



■i6j 



\Qj,Pj) 



(16) 

where Oj — 2ttD Yjl=o Si . They are labeled by the binary 
symbol of the periodic orbit and by the discrete index k 
(k = • • • L— 1). Within the validity of the above approx- 
imations they are orthogonal. The fact that the eigen- 
values are complex reflect the instability of the orbit and 
characterize these states as long lived resonances whose 
approximate width on the unit circle is A. As this width 
is classical (independent of D) these resonances overlap 
significantly with a number of eigenstates \D/2tt which 
is large in the semiclassical limit. 

It is convenient to impose the map symmetries to the 
POM's. The symmetries of the periodic orbits can be 
used to this purpose. The PO of the classical baker map 
can be classified in terms of their invariance under the 
classical symmetries R and T. We characterize this in- 
variance by two integers gr,gt with the value gr = 
or (Jt = for invariant orbits, while ctr = 1 or gt = 1 
if there are two different orbits connected by the respec- 
tive symmetry. As the action S u is invariant under these 
symmetries, the eigenvalues of the POM constructed for 
each S u are degenerate with an associated subspace of 
dimension 2 aT 2 aR . In these subspaces it is possible to 
construct POM's that have the same symmetries as the 
eigenfunctions. Thus central orbits (ctr = a T = 0) with 
v = = V give rise to L states, which automatically 
have the required symmetries. Orbits with either <jr = 0, 
a T = 1 or ctr = 1, gt = give rise to 2L states, while 
non-symmetric orbits gr = 1, gt = 1 give rise to 4L 
states. Some examples illustrating this construction are 



01> 



k \ — 



4) 



001 



1^001011 / 



I<&1> 



I*; 



k,T=±,R=± 
0001011 



(i±r) 

(1±T) (1 ± R) 
V2 V2 



14 



001011/ 



10, 



0001011 



(17) 
(18) 

(19) 

) (20) 



where k = 0, . . . , L — 1. 

Figure [T] shows the Husimi representation of two sym- 
metrized POM corresponding to v = 001 and v = 00101. 
The bottom part shows the distribution of the squared 
overlaps with the eigenstates as a function of the eigen- 
phases. The central dotted line is the Bohr-Sommerfeld 
energy in Eq. [151 It should be clear that this construc- 
tion can be justified for a fixed periodic orbit, and for 
D — > oo because the periodic points of chaotic systems 
are isolated. However if we want these quasimodes as a 
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Figure 2: (color online) top: Husimi representation of 
the scar functions, in the unit square phase space (q,p £ 
[0,1)), over v — 001 with different propagation times 
l$ooi°'™7> (t=0) (left), |^ fc=1 ' te -> (t=2) (center) 

and |$ooi 4 ' fe ~ 1,jR ""> ( t=4 ) ( ri g ht ) for D = 12 0- We plot 
the square of the Husimi function to enhance the structure 
on the manifolds, bottom: Square modulus of the prod- 
uct of these states |$ooi — 1,R ~ ~}\ 2 with eigenstates 
with j = 1, . . . , D ordered by increasing eigenphases. The red 
dashed line represents the central phase Aq^i 



basis for a fixed value of D we need also to consider orbits 
where the assumptions in Eq. [12] are not satisfied. In- 
stead of a diagonalization by means of an explicit Fourier 
transform we have to consider the generalized eigenvalue 

problem det (qi,Pi\B\qj,Pj) - \(qi,Pi\qj-Pj) =0. We 
explore this problem in connection to homoclinic orbits 
in the Appendix A. 



Scar functions 



where S t (s) is the Fourier transform of . When this 
operator acts on a POM it sharpens its quasi-energy 
width while extending the wave packet in phase space 
along the stable and unstable manifolds of the periodic 
orbit. A wave function is thus created that interpo- 
lates between a simply constructed but relatively unsta- 
ble state and a true eigenstate if the propagation time 
is of the order of the Heisenberg time. The interesting 
region is of course when the propagation time is of the or- 
der of the Ehrenfest time. We define the s car functions as 
the short time propagation (oc Ehrenfest time = log 2 D) 
of POM in v with the phase evaluated in e = defined 
at Eq. [G2 



= -P t (2wAZ) |*J 



(23) 



where n is a normalization factor, and the POM is recov- 
ered for t = (|$^ =0}fe ) = Notice that, for unstable 
periodic orbits, the forward and backwards propagation 
of wavepackets on the periodic points lead to significant 
amplitude on the stable and unstable manifolds of that 
orbit. It is then expected [TO, QjJ that these scar func- 
tions will be structures supported on these manifolds and 
having narrower overlaps with definite quasienergy re- 
gions on the unit circle. In this work we fix the phase e 
to the Bohr-Sommerfeld values (A^) and vary the time 
propagation t. 

The Husimi representations of the scar functions over 
v = 001 for different evolution times (t = 0,2,4) , and 
their products with QBM eigenstates are shown in Fig. 
El Note that the coherent states in the Husimi represen- 
tation spread in the stable and unstable manifolds inter- 
fering between each other. Note that the Ehrenfest time 
for D = 120 is approximately ^Ehr ~ 7. 



In the previous section we have seen that the POM 
are quasienergy wavepackets of constant classical width 
A. Narrower wavepackets can be constructed by Fourier 
transfor ming the POM's evolved for a limited time, @,B 
l9l IToL HlLll3j . The resulting states are the scar functions. 

Consider the following operator which depends on the 
time parameter t and the phase e 

00 2 
P t (e) = e~ isl e-^B l (21) 

i— —00 

When the gaussian window is allowed to have infinite 
width (t — > 00), this operator projects on the quasienergy 
eigenstates while for finite t the corresponding width will 
be Ae = 2ir/t In general we have 

D 00 2 

j—l l= — oo 
D 

= Y, 5 ^ £ -^j)^j\ (22) 



IV. NUMERICAL TEST OF THE BASIS 

A. Scar Function Basis 

The QBM can be simplified using the scar functions 
as basis. In fact, it is usefull to choose a set of Ds non- 
orthogonal scar functions, which span the D-dimensional 
Hilbert space, as an overcomplete basis of the QBM. The 
election of the periodic orbits will determine the basis. 
As we want to construct functions on short PO we give 
two different rules to choose the PO of the basis which 
converge to the same basis for long D. 

1) The first choice is to select the PO with shortest 
period which span the Hilbert space. For example, for 
D = 100 we will have Ds — 106 scar functions con- 
structed over PO with period up to L = 6, and for 
D = 140 we have to include PO with period L = 7 and 
therefore D s = 232 (see Table U for D < 472). 

2) The other possible election of the PO could be, for 
a given D, to choose Ds = 2 L G [D,2D) and the short 
PO labeled by binary strings of L digits. 
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dimension 

D s 


L v 


Periodic orbit added 

V 


2 


1 





4 


2 


01 


10 


3 


001 


22 


4 


0001; 0011 


52 


5 


00001; 00011; 00101 


106 


6 


000001; 000011; 000101; 000111; 001011 


232 


7 


0000001; 0000011; 0000101; 0001001; 

nnnm i i . nnm ni i . nnm i ni . nm nm 1 . nm ni ni 
UUUUiii; UUU1U11; UUU11U1; UU1UU11; UU1U1U1 


472 


8 


00000001; 00000011; 00000101; 00001001; 
00000111; 00001011; 00001101; 00010101; 
00011001; 00010011; 00100101; 00001111; 
00010111; 00011011; 00101011; 00101101 



Table I: Periodic Orbits used in basis construction of dimen- 
sion D s to describe the QBM spectrum of dimension D (where 
D S >D). 
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scar functions 



Figure 3: Overlap matrix, K^il^ 



t,k,R,T\\2 



) | , between the eigen- 



states of the QBM ordered by eigenphases and the scar func- 
tions constructed on the PO with L < 6, ordered by increasing 



phase A u for D = 100 (\ipj) on rows and \$J 



*=*Ehr/2, k,R,T 



) on 



columns). The value of the overlap is in grayscale from 
(white) to 1 (black). 



In this work we will choose the shortest period basis, 
but for long D limit both bases, are similar. The im- 
portant fact is that in both cases the maximum period 
involved grows as log 2 D. Fig. [3] shows the overlap matrix 
|(^|$^ ,jR ' T >| 2 of QBM eigenstates (D = 100) ordered 
by eigenphases (l^j), on rows) and the scar function ba- 
sis (D s = 106) (|$£=W2,*,h,t^ Qn columns ) ordered 
by phase value with t = tEhr/2. 

This result, besides restating the fact that scar func- 
tions have narrow overlaps with eigenstates also shows 
that an eigenstate can be pictured as a narrow superpo- 
sition of scars. Each of which is relatively simple to con- 
struct. For example, in the case of D = 102 one of the 
eigenstates of the QBM can be represented by only one 
scar function of period L = 3 with a propagation time of 



Figure 4: Husimi representation in the unit square phase 
space (q,p £ [0,1)) of an eigenstate of the QBM for D = 
102 (left), and its most prominent scar function components 
(right) as given in Eq. [24l 
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Figure 5: (color online) Top: Average participation ra- 
tio ({PR)) as a function of the propagation time for D = 
34, 66, 130 in units of Ehrenfest time; Bottom: VoPRi.b, frac- 
tion of scar functions which have PR < 1.5 in the QBM eigen- 
states basis. 



t = with an accuracy of | (^j l^ooi^ 2 ^ + )l — 

0.87. The approximation of the eigenstate can be im- 
proved adding another scar function of period L = 4. 
Therefore an ansatz state constructed as 



°l 1^001 



)+c 2 |<£ ooi / 
(24) 

with ci c± 0.75 and c 2 — 0.48e 20117 has an accuracy of 
| (ijjj | <E> ansatz) | — 0.92. Figure [4] shows the husimi func- 
tions of the eigenstates and both scar functions consid- 
ered in this example. 

The scar functions have superposition with a small 
number of eigenstates. This can be quantified in terms 
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D 

Figure 6: (Color online) Average participation ratio ((PR)) 
as a function of the dimension of the QBM for the POM and 
the scar functions with t — tEhr/2 and t — tEhr- 



of the average participation ratio defined as 



(25) 



L=l \ 3=1 



where (PR) G \l,D). Notice that the random matrix 
theory prediction for a generic basis is (PR) = D/2. 

The (PR) decreases with the propagation time and 
converges to 1 in the Heisenberg time limit as was ex- 
pected. However, the interesting region is for times 
of the order of the Ehrenfest time. In Fig. [5] (top) 
we show (PR) as a function of propagation time for 
D = 34, 66, 130. In Fig. [5] (bottom) we also show the 
fraction of scar functions which have PR less than 1.5 
as a function of time. Notice that for D = 130 and at 
t ~ tEhr the average number of eigenstates in a scar is 
about 3 and 20% of the scar states have a < PR > less 
than 1.5 meaning that they are almost pure eigenstates. 

The average participation ratio as a function of the di- 
mension of the QBM is shown in Fig. [6] for the POM and 
the scar functions with t = tEhr/2 and t = tEhr- For the 
limited range of values available the (PR) seems to grow 
linearly with D but with a slope significantly smaller than 
the D/2 random matrix value. The reduction is similar 
(and more important) than that obtained in [l|, |2(]. 



V. SCAR APPROXIMATION BY HOMOCLINIC 
PERIODIC ORBIT MODES 

The construction of the POM basis is a simple an- 
alytical formula only involving as classical input short 
periodic orbits. The scar function requires in addition 
the forward and backwards propagation of the POM. In 
this paper for simplicity this propagation was obtained 
exactly by a matrix multiplication. It should be clear, 
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Figure 7: (Color online) Intensity I„ = |($* E W2. fc - 
between the scar function over u = 01 with t = tEhr/2, and 
the homoclinic POM over \x = (01) s with s = 3, 5, 7. k v = 
is shown on top and k u = on bottom, where respective 
phases are drawn with solid black lines. The dimension of the 
Hilbert space is D — 128 and the products are ordered by 
increasing phase (0) of \i for the 2s + 1 values of k^. 



Figure 8: Husimi representation in the unit square phase 

space of the scar function |3v Ehr/2 ' ) on the PO v = 01 
(left). On the right, the Husimi representation of the POM 
|^o,fc=i2^ cons t mc ted on the homoclinic orbit fi = (01) s 
(with s = 7) which approximates the scar function. We plot 
the square of the Husimi function to enhance the structure 
on the manifolds. Those states are the scar function and the 
POM approximation of the scar of Fig. [7] respectively for 
D 128. 



however, that accurate semiclassical expressions for this 
propagation could be available as long as the times in- 
volved remain bounded by the Ehrenfest time. In this 
section we illustrate a different approach that replaces the 
need for this propagation by the construction of POM's 
on long periodic orbits homoclinic to short ones. 

The homoclinic POM are simply the POM construe- 
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tion described in Section IIII Al over PO of the form 
\i — (v) s h with period = sL u + L^. This PO re- 
sembles the homoclinic trajectories of v with the excur- 
sion h represented by the infinite string . . . vwhvvv . . ., 
and converges to it in the long s limit. These orbits for 
large s accumulate near stable and unstable manifolds of 
v and therefore can mimic the building up of amplitude 
produced by the propagation. 

A long POM with period will generate states 

k 

with phases A£ uniformly spaced on the unit circle. The 

k h 

state with the phase closer to A* u of the short v or- 
bit will have the largest overlap with the scar state. An 
example is shown in Fig. [71 for the product between the 
scar function on v = 01 for both values of k u (k = 
on top and k = 1 on bottom) with t = tEhr and the 
homoclinic POM over fj, s = (y) s h for all values of k^ 
(I u = K^Ehr,^ |$^7°' fe/x )|) as a function of the phase 
which take the discrete values A* . In this case we choose 
the shortest homoclinic excursion h = and values of 
s = 3, 5, 7. In this example acceptable approximations 
(I v > 0.9) are reached for values of s no longer than 7. 
In Fig. [8] we compare the husimi representations of the 
scar function on v = 01 with k = 1 and t = ^Ehr/2 an d 
its best approximation with the homoclinic POM over 
fi s = (01) s with s = 7. Note that the homoclinic POM 
structure looks similar to the evolution in time of the scar 
since it spreads the stable and unstable manifolds of the 
PO. 



VI. CONCLUSIONS 

We have shown that a set of states constructed on the 
periodic orbits of the QBM provides a way of describ- 
ing its eigenfunctions which significantly improves (in 
terms of the average participation ratio) the description 
in terms of a generic basis. For the QBM this improve- 
ment is similar to that obtained by Lakshminarayan us- 
ing the Hadamard basis [l| or by the essential baker map 
0| . However the method does not use any special struc- 
ture of the chaotic map and should be generally applica- 
ble to other chaotic maps [21] as long as the properties 
of a limited set of periodic orbits are known. 
The construction is still far from a full semiclassical anal- 
ysis. The basis that we constructed is non-orthogonal 
and requires to consider each primitive orbit as non in- 
teracting with each other periodic orbit. The semiclas- 
sical calculation of this interaction would be needed to 
calculate the amplitudes that describe the eigenstates in 
terms of scar states. Calculation of this type has been 
performed for the billiard stadium [TO, Hlj l2p , an hyper- 
bolic Hamiltonian [23], and the cat map [211 ]. 

It is also important to mention that the aim here is 
to reconstruct the unitary dynamics of the map in terms 
of pure states constructed as interpretations of gaussian 
packets with complex coefficients derived from the clas- 
sical action. A different construction based on incoher- 



ent superposition of densities placed on periodic points 
would allow a similar reconstruction in terms of the Liou- 
ville dynamics. We postpone this aspect of the problem 
for future investigations. 

We thank Eduardo Vergini for interesting discussions. 
Partial support by ANPCyT and PIP6137 of CONICET 
are gratefully acknowledged. 



APPENDIX A: POM Diagonalization 
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Figure 9: (color online) Top: Real and imaginary part of 
the solutions of the generalized eigenvalue problem for the 
PO v — 01 and the homoclinic orbit family /x — (01) s for 
s = 3, 4, 5, 6, 7. with D = 50. These solutions converge to the 
exact eigenvalue which lives in the unit circle (in solid line). 
Bottom: Husimi functions in the unit square phase space of 
the states associated with the eigenvalues marked with • (/j, = 
01), + (/x = (01) 7 0) and x (the exact eigenstate of the map) 
on the left, center and right respectively. 



The periodic orbit modes (POM) constructed in Sec- 
tion [IIlAl can also be seen in a different light. Consider 
the generalized eigenvalue problem 



det 



((7*l%i>-^7i|7i>) = 







(A-1) 



where are coherent states placed on M periodic 

points of the map. The operator Y^f l7i)(7i| can be 
considered as an approximation to the resolution of unity 
in the D-dimensional Hilbert space, characteristic of co- 
herent states as long as M is large enough (> D) and 
the points cover uniformly the phase space. This latter 
condition is satisfied for chaotic maps on account of the 
Hannay-Ozorio de Almeida sum rule [24|. Under this 
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conditions Eq. IA-11 provides D exact eigenvalues and 
M — D zero eigenvalues [25| . 

POM are obtained with the two assumptions: 

a) Different primitive orbits do not interact. 

b) Eq. [12] is satisfied. 

With these conditions the eigenvalue problem is solved 
by a simple Fourier transform involving the different 
points on an orbit. The eigenfunctions are the POM's 
of Eq. [16] and the eigenvalues are complex and given by 
Eq. n~Q 

If we retain only retain the assumption that periodic orbit 
do not interact we can refine the construction of POM by 
considering the limited diagonalization (Eq. lA-lj) where 
now \ are the points on a primitive orbit (and its image 
under R and T if necessary). This refinement is impor- 
tant for long orbits of the type fj, = v s h and fj, = v^hv^ 
which mimic homoclinic and heteroclinic orbits and ac- 
cumulate their points close to short primitive orbits. For 



these orbits it is not correct to assume the coherent states 
as being even approximately orthogonal and therefore the 
generalized diagonalization is necessary. The eigenvalues 
divide in two sets. Some converge rapidly to very small 
values, while there are always some of them which con- 
verge to the unit circle. 

The advantage of this construction is that the "good" 
eigenvalues are closer to the unit circle and therefore 
represent modes that have a longer lifetime. Moreover 
spurious eigenvalues that are produced by the superpo- 
sition of many almost equal coherent states are rapidly 
eliminated. The disadvantage is of course that the con- 
struction is not analytic and requires a diagonalization. 
We give an example for the family of orbits \i s = v s h. 
We have diagonalized Eq. IA-11 for periodic points on the 
orbit i± s (s = 3,4,5,6,7) and compare with the exact 
eigenstate of the QBM for D = 50. 
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